Optimal State Discrimination Using Particle Statistics 



S. Bose 1 ' 2 , A. Ekert 3 , Y. Omar 4 ' 5 , N. Paunkovic 4 and V. Vedral 6 
1 Institute for Quantum Information, California Institute of Technology, CA 91125 
2 Department of Physics and Astronomy, University College London, Cower Street, London WC1E 6BT, United Kingdom 
3 Centre for Quantum Computation, DAMTP, University of Cambridge, Cambridge CB3 OWA, United Kingdom 
4 Centre for Quantum Computation, Clarendon Laboratory, University of Oxford, 
Parks Road, Oxford 0X1 3PU, United Kingdom 
5 Centra de Fisica de Plasmas, Instituto Superior Tecnico, P-1049-001 Lisbon, Portugal 
6 Optics Section, The Blackett Laboratory, Imperial College, London SW7 2BZ, United Kingdom 

(September 8, 2003) 



We present an application of particle statistics to the prob- 
lem of optimal ambiguous discrimination of quantum states. 
The states to be discriminated are encoded in the internal de- 
grees of freedom of identical particles, and we use the bunch- 
ing and antibunching of the external degrees of freedom to 
discriminate between various internal states. We show that 
we can achieve the optimal single-shot discrimination proba- 
bility using only the effects of particle statistics. We discuss 
interesting applications of our method to detecting entangle- 
ment and purifying mixed states. Our scheme can easily be 
implemented with the current technology. 



Recently, there has been an emerging interest in the 
use of particle statistics (both bosonic and fcrmionic) for 
quantum information processing [1,2]. In fact, it was 
shown that a useful task such as entanglement concen- 
tration could be accomplished, even if non-optimally, us- 
ing only the effects of quantum statistics, without the 
need for any other interactions [2]. The above inves- 
tigations differ significantly from some previous sugges- 
tions, where either anyonic statistics [3] or the effects 
of electronic statistics in conjunction with other interac- 
tions [4,5] were used for quantum information processing. 
Schemes using only particle statistics [1,2] would be very 
useful for tasks implemented with identical particles that 
interact very weakly or not at all with each other, such as 
photons or neutrons. This weak interaction can be bene- 
ficial for information processing as it may reduce the un- 
wanted coupling to the environment. Such schemes are 
also extremely general in the sense of being independent 
of the actual particle species. It is however not known 
whether such schemes can accomplish quantum informa- 
tion processing efficiently. Here we present a particu- 
lar quantum information processing task involving two 
qubits and show that it can be performed optimally using 
only quantum statistics. Moreover, we point out how the 
task of discriminating quantum states can be applied to 
detecting entanglement and purifying mixed states. We 
also discuss how to generalize these tasks to N qubits 
and argue that quantum statistics could be used to per- 
form even this generalized task optimally. While the two 



qubits and other small iV qubits versions of our protocol 
can be tested with photons, electrons, neutrons or atoms, 
the large number of qubits versions could have interesting 
implementations in optical lattices [6,7]. 

One of the striking aspects of quantum mechanics is 
that it is not possible to perfectly discriminate between 
two states unless they are orthogonal. Suppose someone 
prepares two qubits encoded in the internal degrees of 
freedom of two identical particles — say, in the spin of 
two electrons or the polarization of two photons — in one 
of the following two possible states: 

• spins aligned (parallel) and pointing in an arbitrary 
direction; 

• spins anti-aligned and pointing in an arbitrary di- 
rection. 

We will assume that the two states are equally likely for 
the sake of simplicity, but all the presented results will 
be valid for any a priori distribution. Note also that, 
unless otherwise stated, whenever we mention spin we 
will actually refer to any two-dimensional internal degree 
of freedom, be it for fermions or for bosons. As we do not 
have any knowledge of the direction of alignment in both 
of the above cases, the overall states are mixed. They are 
described respectively by the following density operators: 

P 2 = -^[(m\n){n\®\n){n\, (i) 

and 

er 2 = — / dfl\Q){Q\ <g> |^- L )(f7- L |. (2) 
4tt J 

Here, the subscripts indicate that we are considering two 
particles. The ket |0) represent the spin- up state along 
the axis defined by the angle f2, while \£l ± ) is the or- 
thogonal state spin-down. Thus, p2 represents an equal 
mixture of spins aligned along an arbitrary axis in space, 
while 02 is the equal mixture of anti-aligned spins of two 
spin-1 particles. It is impossible to discriminate between 
these states perfectly because they are not orthogonal. 

Optimal results are known for the discrimination of 
any two given quantum states r\ and r)' [8] . The maximal 
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probability of ambiguously discriminating between two 
a priori equally likely quantum states in a single-shot 
measurement is given by the Helstrom formula: 

PH(v,v') = l + \Tr\r 1 -r 1 '\. (3) 

We now present a procedure for discriminating be- 
tween p2 and (72, both for fcrmions and bosons, based 
only on the effects of particle statistics. To use these ef- 
fects we interfere particles at a beam splitter. Here we 
use beam splitter in a generic sense, referring not only 
to the common optical element (partially silvered mir- 
ror) used with photons, but also to any device presenting 
an analog behavior for other kinds of particles, as al- 
ready suggested for electrons [9]. In order to distinguish 
the states we rely on path measurements that discrimi- 
nate between bunching and antibunching, and that are 
performed on our particles after letting them pass simul- 
taneously through a 50/50 beam splitter. Note that in 
such balanced beam splitters two indistinguishable parti- 
cles will always bunch if they are bosons [10], and always 
antibunch if they are fermions [9] (sec also [11]). 

For fermions, our guess in the case of the antibunching 
result is that the input state was p 2 , while in the case of 
bunching is that it was u 2 . The probability of success of 
our procedure is then: 
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Pbs(P2,t 2 ) = 4. (4) 

This probability can easily be calculated by noticing that 
the only case for which our guess could be incorrect is 
when we have the antibunching result. Then, we con- 
clude that the input state was p 2 , while it could have 
actually been 02. On the other hand, in the case of the 
bunching result we know for sure that the input state was 
(T 2 , since — according to the Pauli exclusion principle — 
two particles with aligned spins cannot end up in the 
same output arm of the beam splitter. When our input 
state is (T2, the antibunching happens with probability 
1/2, giving in total a probability of incorrect inference of 
1/4. In the case of bosons, our protocol is exactly the 
opposite of the fermionic one, but yields precisely the 
same efficiency: this time the antibunching results stand 
for cr 2 , whereas the bunching ones stand for p 2 , but the 
probability of success coincides with Eq. (4). 

Interestingly, the Helstrom formula gives the same re- 
sult for the maximal probability of discriminating be- 
tween these two states: 

P H {P2,G2) = \. (5) 

We can thus conclude that our procedure is optimal for 
both fermions and bosons. 

Helstrom's probability can still be achieved in other 
cases using the effects of quantum statistics. We now 
introduce a case of special interest that can be applied 



to other tasks such as entanglement detection and state 
purification, as will be shown later in the paper. Suppose 
that we have to discriminate between the following two 
states: 

• spins aligned (parallel) and pointing in an arbitrary 
direction (the same as in Eq. (1)); 

• each spin in the maximally mixed state. 

The latter state is represented by the following opera- 
tor: 

t 2 = \l® 2 = \ (|0)(0| + |1)(1|) ® (|0)(0| + |1)(1|) , (6) 

where |0) and |1) are any two orthogonal spin states. 

The strategy now is exactly the same as before. It re- 
lies on the fact that if the state is p 2 the particles can, 
due to their indistinguishability, give only one result (an- 
tibunching in the case of fermions or bunching in the case 
of bosons), while if the input state is T2 both results are 
possible. This time, the probability of success is: 

Pbs(P2,t 2 ) = 1, (7) 

and this coincides with the Helstrom result 

Ph( P 2,t 2 ) = 1, (8) 

so our procedure is optimal in this case too. The calcu- 
lations leading to the above results are analogous to the 
ones in the previous discrimination case. 

One of the most interesting applications of our ap- 
proach can be found in the detection of entanglement. 
We illustrate this in the case of pure states of two parti- 
cles. Suppose, for example, that we are to discriminate 
any maximally entangled state of two qubits from any 
disentangled (product) state. As before, the qubits are 
supported by the internal degrees of freedom of two iden- 
tical particles (labelled, say, A and B). Suppose in addi- 
tion that we are given two identical copies of the state. 
In order to detect entanglement, we take the same par- 
ticle (either A or B) from each pair and interfere them 
at a beam splitter, as shown in Fig. 1. The crux of the 
argument is that if the state is entangled, then the re- 
duced states of these particles will be maximally mixed 
as in the state r 2 . On the other hand, if the state is 
separable, then the two interfering particles are in the 
state p 2 . This is the same as in our discrimination pro- 
cedure above, and so there is a probability of 5/8 to de- 
tect entanglement. This example can be generalized to 
other entangled pure states, and, more interestingly, to 
more particles. We note that there is a close analogy be- 
tween this method and our entanglement concentration 
scheme in [2]. In particular, if the two states were less 
than maximally entangled, then by detecting entangle- 
ment we would actually also amplify it (see [2] for more 
details). 



2 



Another interesting application is in mixed state purifi- 
cation, as in [12]. Suppose that we start with two qubits, 
each in some mixed state. We would like to make the 
state of these qubits purer (in the sense of having lower 
linear entropy), but also to preserve their original direc- 
tion in the Bloch sphere. The optimal way of doing so (as 
proven in [12]) is to project the joint state onto the sym- 
metric subspace, in which case the resulting mixed state 
is purer and yet preserves the original direction. If the 
projection is unsuccessful, the qubits are thrown away. 
This is exactly the same as our probabilistic discrimina- 
tion with a beam splitter. 

We would now like to investigate the generalization 
of the above results to N particles. For this, we use a 
generalized (N-port) balanced beam splitter, as shown in 
Fig. 2, which acts only on the spatial degrees of freedom 
of the input particles. This action is given by a unitary 
matrix Un, with elements: 
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i^f(m-l)(n-l) 



(9) 



(Note that there exist alternative descriptions of bal- 
anced multiport beam splitters [13].) The square of the 
norm of each element in the matrix represents the proba- 
bility that the particle in the m-th input arm of the beam 
splitter ends up in the n-th output arm. Since all these 
elements have norm -^=, we have a representation of a 
balanced N-port beam splitter. 

The aim is now to discriminate between the N-particle 
generalizations of pi and ti. Those states are given by: 



pn = ^-J dn(\n){n\f N 



and 



™=2AT(|0><0| + |1><1|) 



(10) 



(11) 



To calculate the Helstrom probability, one has to diago- 
nalize the matrix pn — tn- This turns out to be straight- 
forward once we notice that pn can also be represented as 
an equal mixture of all possible symmetric states within 
the basis {\S l } : i = 1, . . . , N + 1 = \S\} of N qubits. 
Then, we have: 



Pn 



\s\ 

' ' i— 1 



(12) 



It is now easy to calculate the Helstrom formula by ex- 
panding pn in a basis consisting of the union of a basis 
of the symmetric sub-space and a basis of its orthogonal 
complement. The result is: 



Ph(pn,t n ) = 1 



(N + l) 



(13) 



Alternatively, we can calculate the average probability of 
success to distinguish states pn and tn using the follow- 
ing expression: 



Ph(pn,t n ) = il + ip. 



(14) 



Here, the 1/2 factors refer to the fact that the two states 
are prepared with equal probability. The term 1 comes 
from the fact that the state pn, supported on the sym- 
metric subspace (of dimension d$ — N+l), is always 
identified reliably as such, and p is the probability of 
identifying the other state, tn- Since tn is maximally 
mixed, it is uniformly distributed over the whole space of 
N qubits (of dimension d — 2 N ). In this case p = 1 ~ d s , 
which after substitution in Eq. (14) gives Eq. (13) right 
away. 

The problem to apply our discrimination scheme to N 
particles is that it becomes exponentially hard to cal- 
culate Pbs{pn ,t~n) as N increases. Moreover, it is not 
clear which inference strategy should be followed. For 
fermions, the natural generalization seems to be to as- 
sociate the antibunching results with pn and the others 
with tn- For bosons, on the other hand, a more subtle 
strategy may be needed. This is because, loosely speak- 
ing, for bosons there is no clear analogue of the Pauli 
exclusion principle. Furthermore, even without having 
the complete calculations for N > 2, we would like to 
emphasize the remarkable fact that the Helstrom prob- 
ability Ph{pn ,t~n) is equal to the probability of success 
of a fermionic beam splitter strategy described above, if 
calculated under the assumption that the particles are 
classical (i.e., always distinguishable by some arbitrary 
label), but obey a constraint equivalent to the Pauli ex- 
clusion principle (not allowing more than two particles 
in the same internal state to share the same output arm 
of the beam splitter) . The overall probability is then cal- 
culated by summing up the probabilities of all the possi- 
ble outcomes rather than the amplitudes, as it would be 
done in the quantum case. However, in the case of three 
fermions (N = 3), we have performed the full quantum 
calculations (i.e., taking properly into account the effects 
of statistics) for a three-port balanced beam splitter and 
obtained: 



Pbs(P3,t 3 ) = -, 



(15) 



which is equal to Ph(P3, T3). We believe this result of ob- 
taining the optimal discrimination probability using the 
effects of particle statistics (in multiports) can be gener- 
alized to an arbitrary N, both for fermions and bosons, 
and we continue research in this direction. For now this 
remains a conjecture. In an optical lattice, with one par- 
ticle in each lattice site, a multiport beam splitter could 
probably be simulated by dissolving N potential wells and 
then creating a new set of N wells [7] . Of course, if the 
particles interact, then the effective beam splitter will be 
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modified, and here we only point out the plausibility of 
creating multiport beamsplitters (or multiparticle inter- 
ference) in an optical lattice. 

In this paper we have shown that it is possible to per- 
form an optimal quantum information processing task 
using only the effects of particle statistics. In particu- 
lar, we have presented a strategy for discriminating be- 
tween two non-orthogonal states of two qubits (encoded 
in the internal degrees of freedom of identical particles) 
using beam splitters. We have considered two discrimi- 
nation scenarios and in each of them our strategy differs 
(symmetrically) between fermions and bosons, but offers 
the same efficiency. We also pointed out how our dis- 
crimination scheme can be applied to detecting entangle- 
ment and purifying mixed states. In addition, we have 
calculated the Hclstrom probability for N qubits in one 
of our discrimination scenarios. We have shown by ex- 
plicit calculation that this probability can be achieved in 
a fcrmionic three-port beam splitter strategy and that it 
is the same as the fermionic strategy for general N if the 
fermions are considered as classical particles that obey 
the Pauli exclusion principle as the only additional con- 
straint. An advantage of our method is that it can also 
be easily implemented with the current technology. 

Our work suggests a number of interesting research di- 
rections. One problem is to prove the optimality of the 
beam splitter strategy in the case of N qubits and its ap- 
plication to multiparty entanglement detection. This, we 
hope, will answer the question of weather the symmetry 
between fermions and bosons in our strategy will be pre- 
served for a generalized beam splitter. It may also lead to 
a simple and physically intuitive selection principle gov- 
erning bosonic behavior. Another possible direction is to 
classify all the pairs of states that can be optimally dis- 
criminated with our scheme. Finally, our results suggest 
that it would be worth to further explore the role of parti- 
cle statistics in quantum information tasks, in particular 
in efficient quantum computation. 

After the conclusion of this paper, an interesting re- 
lated work was pointed out to us [14]. The authors would 
like to thank G. Castagnoli on many useful discussions 
on the link between computational complexity and par- 
ticle statistics. His ideas have stimulated much of the 
work in this paper. We also would like to thank F. A. 
Bovino and P. Varisco for useful discussions on the role 
of particle statistics in quantum information processing. 
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FIG. 1. This figure represents the set-up for our entangle- 
ment detection scheme. We have two equal pairs of identical 
particles and consider their internal degrees of freedom. Both 
pairs are in the same pure state: either separable or maxi- 
mally entangled. We take the same particle (for instance, B) 
from each pair and interfere them at a 50/50 beam splitter. 
If the states are disentangled then the particles are indistin- 
guishable and, depending on the statistics, will either only 
bunch or only antibunch. Otherwise, if the states are maxi- 
mally entangled, the particles are maximally mixed, meaning 
that they can be (probabilistically) distinguished and hence 
the statistics does not influence their behavior. This entan- 
glement detection is a particular instance of our state discrim- 
ination scheme discussed in the paper. 
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FIG. 2. This diagram represents a multiport beam splitter 
with N inputs and N outputs. The overall output state de- 
pends not only on the input, but also on the statistics (either 
fermionic or bosonic) of the identical particles involved. We 
have labelled two arbitrary ports, the m-th input port and 
the n-th output port. 
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